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ABSTRACT 


A computer program for Simulating submarine motion in 
Six degrees of freedom is developed. The Simulated 
Submarine iS given a capability of shifting ballast. An 
automatic pitch and depth control is designed for the 
Submarine simulated using optinal control theory. With the 
depth and pitch in automatic control a trim error Signai is 
developed by comparing the parameters of the plant to that 
of a linearized model. This error is used to implement an 
automatic trim control that shifts ballast within the 


Submarine to achieve a neutcal trin. 
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I. SUBMARINE CONTROLS 

| The submarine control problem differs from that of the 
surface ship in that the submarine translation in three axes 
and rotation about each of these axes must be considered. A 
typical submarine implements control about these axes with 
the use of movable plane control surfaces, shifting of 


ballast water and the propulsion systen.| 


The submarine operates and is controlled in all six 
degrees of freedom. To maneuver, three sets of plane 
Surfaces, the propuision system, anda means of shifting 
ballast are used.) This thesis considered only control of 
the plane surfaces andthe shifting of baliasSt water. In 
addition , of the three plane surfaces available only the 
fairwater planes and stern planes were automatically 
@ontrolled, 1.é€. no automatic rudder. The three modes of 
control were depth, pitch and bailast controi. Course 


control was not considered. | 


| Depth control consists of the sensing of any error 
between ordered depth and pitch and the actual depth and 
pitch. The planeS were used to implement the control and 
bring the submarine to the ordered depth and pitch thereby 


reducing the error to zero or near zero. 


Ballast control is the Sensing of other than neutral 
trin ova the submarine. Because or surrounding water 
temperature changes, taking sea water in the ship or pumping 
it out, the submarine may not weigh the same as the volune 
of water it displaces. Also too much water forward or aft 
in the ship initiates a turning moment around the Y axis 
that will eventually have to be corrected by the planes. 
This condition must first be detected so a corrective action 
can be applied by shifting ballast within variable ballast 


tanks located inside the Submarine. 









The purpose of controiling depth is obvious. pitch 
control 1s required largely because it is coupied to depth 
control and it is sometimes desired to maintain pitch at 
Other tnan a zero angle. A means of ballast control is 
required to surface the ship, maintain depth under certain 
conditions, allow depth control to be accomplished within 
the limits of the planes' capabilities and to reduce drag 
caused by a steady state non zero planes angle. fhe 
corrective action for the last two out of trim conditions is 
normally accomplished by use of a pump to shift water into 


Or Out OL variable ballast tanks in the submarine. 


The method of achieving an automatic control was to 
first develop an automatic controller that maintained an 
ordered depth and pitch. The controller was to counteract 
external forces by movement of the fairwater and stern 
planes. Tnese forces can be caused by an out of trin 
condition or by the ship dynamics which at difterent speeds 
would cause the submarine to change depth in varying degrees 
plus such things as sea state, etc. The controller designed 
is unique to the submarine considered in ref(1) but the 
method is general and may be applied to most situations. 
With the depth controller in operation the automatic control 
of ballast was then solved. Again the solution 1S unique to 
the Submarine considered. Finally with all automatic 
controls in operation various runs were made to determine 


the capability of the systen. 


However, betore a controller deSign can be started a 
computer model one the submarine aust be developed. 
Reference 1 containS a FORTRAN program used by NSKDC. A 
program ianguage that waS more easily nanipulated in terms 
of filters, integration routines, etc. was DSL. To use DSL 


the equations of motion would nave to be solved. 
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Il. MODELING A SUBMARINE IN A DIGITAL COMPUTER 


——— —— eee ee oe oe — a co eee ee eee ee ee =—— ame = Reo ame So eee ee rs ce eee ee 


A. EQUATIONS OF MOTION FOR A SUBMARINE 


The equations of motion in six degrees of freedom are 
derived by eguating the summation of the forces and moments 
to zero. 

2 Forces=0 
2 Moments=0 
Where the vector force 
F = d/fdt Mfomentun 
=d/fdt (mVe) = md/dt Ve 
m is the mass and Ve is velocity 
Ve = 1U+jV+kW 


mae direction of positive motion is indicated in fig. II-1. 


Moment d/fdt (Angular momentum) at the center of gravity. 


Moment = d/dt (1PIxxt+jolyy+kRizz) 


Where Ixx is the moment of inertia about the xX axis etc. 
The direction of positive motion is shown in Fig. II-2. 
A simplified derivation cf the equations of motion is found 


mm ref. 2. 


The eguations are referred to the body, i.e. the origin 
of the ccordinate system is located in the body and at the 
center of gravity. Figure II-3 is a submarine placed on the 
coordinate axiS with arrows pointing in the direction of 
poSitive motion. In addition the movement or the rudder, 
Stern planes and fairwater planes are also depicted on Fig. 
II-3 witn the arrows shown pointing in the direction of 
positive motion. The specific equations or motion for a 
Submarine were developed for NSRDC and presented in ref. 3. 


If the coordinate origin 1S put at the center of gravity 


11 





these eyuations become the the equations of ref. 1. The 
equations of motion for a submarine in six degrees of 
freedom are repeated from ref. 71 in the Appendix. For use 
in the ccmputer programs that follow the equations 


containing I ,I and I had both sides of the eguation 
x y Z 


S 
divided by 1, Where 1 1s the ship Length. In the 


3 
equations containing m, both sides were divided by 1.f 


waS assigned the value of 2. 


The eyuations contain forcing tecms for the planes and 
rudder, shifting of ballast water and the propeller. The 
propeller input 1S a function fitted to a curve. lee 
contains the command speed and actual speed as ratios. 
Normally there are different sets of coefficients for 
different propelier modes such as packing, slowing down etc. 
The only mode considered was forward propulsion at a 
constant speed. Tt was _noted— thatthe actual speed-was 
always sligntly lower than the ordered ~speed. Uc (commana 


speed) and U(actual speed) are in teet per second. 


All angles are in radians. Mass is in slugs and 


@istance in feet. 


B. COMPUTER SOLUTION TO THE 8QUATIONS OF NOTION 


To solve the reSuitant differential equations ref. 1 
partitions the equations into different sums. The equation 
for U& contains no other acceleration term so it is solved 
separately. V,P and R eyguations all contain mutual 
accelerations so they are solved using Cramers Rule for a 3 
x 3 system. Wand Q contain mutual accelerations so they 


are grouped and solved as a 2 x 2 system. 


Z 





The approach chosen in this thesis was to solve the 
whole system of equations at once using Cramers Rule for a 6 
x 6 system. This was more methodical and the determinants 
were solved only once in the beginning of the program using 
the "INITIAL" capability of DSL. 


The general form of the equations 1S, using two degrees 


of freedom for an example 


Wn 
Ub 
T 
Uhy 
i 
0X 


1 1 12 2 2 
< ete) | 2 2 => iil 


Mh Sth M S+N 
2 22 


K K 
12 21 
2 
xX = I 
S wv nw 
K K : 
12 22 
defining 
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I K S§S 
1 21 
2 
I K §S 
2 22 
es 
1 
4 
S &B 
2 
Kk 5S ie 
11 1 
2 
nS I 
12 2 
= 
2 
& 
o> 4S 
Noting the cofactor of K a = ae and the cofactor of K is = 
ca 
-K etc. 
21 
*m = i cofactor K + I cofactor ik 
1 L Li 2 12 
2 
Say as 
m™ = ~£ cotactor Kk + {| cofactor K 
2 2 22 1 21 
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The same general approach can be used for higher order 
systems with the number of cofactors taken equal to the 
number of variables to be solved. In the eyuations of 
motion I will be a number of time varying nonlinear 


functions. 


Reference 4 used the library subroutine DTERM in 
modified form to solve similar eguations for the CSMP 
language. This method was further modified to work with DSL. 
Six by six determinants were solved by the program to get a 
solution for the equations in the forn of the previous 
paragraph. By substituting columns and rows of ones and 
zeros in the row aud column where the variables I should be 


the cofactor was found. For example in a3 by 3 case 


¢) 
0 K K = K kK - K_ kK = cofactor K 
22 32 22°34 Jg2 23 11 
0 K K 
23 ye! 
K 0 K 
11 31 
K 0 K =m Kk - K kK = cofactor K 
12 32 le 31 114 32 23 
0 1 0 


Program #1 was used to solve for the cofactors in this 
manner. Once the cotactors were found subseguent prograns 
had the cofactors read in aS constants. In the program the 
cofactors are constants and the variable I 1s dependent on 


time and the dynamics of the submarine. 
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C. SUBMARINE SIMULATION 


The computer used for the Simulation was the IBM 360/67 
located at the W.R. cChurca Computer Facility, Naval 
Postgraduate Scnool. The programming language was the Naval 
Postgraduate School version of DSL/360. Computer progran #1 
is a Simulation of a submarine in motion at 15 knots. The 
coefficient values used are obtained in ref. 1. AS was 
previously discussed the eguations of motion are derived in 
the body fixed system. A trigonometric conversion of the 
velocities iS used in the auxiliary equations (appendix) to 
transform the equations to the eartn reference system with 


mie Origin at x(to), y(to) and z(to). 


Reference 17 contains several maneuvers used to check 


out the eguations and coefficients. One wtaneuver is the 

vertical plane overshoot where the stern planes are 
0) 0 

deflected to +20 ; when the ship is pitched to -5 the 


stern piane angle is reversed. Program #1 executes the sane 
Bomeuver with fig. £fi-4 and II-5 showing depth,pitch and 
stern plane angle. The results compared favorably with those 


me ref. 1. 


The original equations contained forcing terms that 
were to represent water being blown from the main ballast 
tanks. In the programs developed for this thesis the 
forcing terms represented water being moved in variable 
ballast tanks with the capability of shifting ballast 


between tanks or between any tank and sea . 


To control the trim of a submarine a means of shifting 
ballast water must be provided. For the computer Simulation 
a tcrigd pump must also be Simulated. Also trim tanks must be 
positioned in the Ship model. The position of the tanks was 


Mie daloitrary Choice wWith Some traditional basis. The tanxs 
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were poSitioned as indicated in fig. II~6. 


Many options are available in the shifting of variable 
ballast water. The options depend on the number of tanks 
and their individual purpose. To simplify the computer 
Simulation the options wiil be decreased considerably. Only 
three tanks will be considered. The auxiliary tank is placed 
at the center of gravity and will be used to adjust any 
Overall mass error. That is a weight differential between 
the snip and the water it displaces. The forward and after 
trim tanks wiil be used only to adjust ballast within the 
Ship. The result will allow negative vaiues for the content 


of a tank. Tae options are, 


1. Heavy- pump auxiliary tank to sea. 

Peeeieght—- £1ood auxiliary tank from sea. 

3. Heavy forward- pump from forward trim tank to after 
mrim tank. 

i tere OrwaLra- pulp from after trim tank to forward 
eam tank. 


Any ballast correction will be accomplished by first 
correcting the overall mass error then achieving a forward 


and aft trim. 


Using the pump fron Program #2 with Program #1 the 
effects of shifting ballast was demonstratea in the 
following simulations. The planes are put at zero and not 
movea during the run. Figure II-7 is the resultant pitch 
change when 2000 pounds of ballast is flooded into the 


auxiliary tank at a@ rate ot 3860 lbs./fmin. Because of the 


0 
rap dynamics a 2.7 up angle is taken even though the 
mixtiiary tank is located at the center of gravity. With 


the planes left at zero the pitch angle remains constant at 
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his angle. (Figure Ii-8)produces a surprising result, with 
2000 pounds ~flooded in the ship went up 240 feet in 3.5 
minutes. (A negative number for depth indicates up and a 
positive number indicates down.) However keeping in mind 
that the submarine assumed an up angle the resultant depth 
change makes more sense. Figure II-9 was plotted with > 


ballast pumped out of auxiliary tank to. Sea. The ship 


0 
assumes a .9 degree up angle. The depth is seen in Fig. 


IZ-10 to change by 135 feet in 3.5 minutes. Prgure —il= 1 
helps to resolve the apparent problem. In fig. II-11 no 
ballast was transferred yet the ship assumed an up angle of 
Significant degree with the planes on zero. Reference 5 
explains that most military Submarines will generate a 
positive pitch due to the ship structure if the planes are 
zeroed. The submarine in ref. 1 ais fictitious and the 
coefficients were randomly selected put it still appears to 
have some of the expected characteristics. ATELY ara sklore 
II-12 and 13 were generated with ballast shifted from after 


beam to fcrward trin. 
The model performed as expected and the trim pump and 


planes produced the desired result. fhe next problem solved 


was the designing of automatic controis. 
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Figure II-1 


Axes Definitions 





Figure II-2 


Axes Definitions 
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Curve 1, Stern Plane Angle vs. Time 
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III. AUTOMATIC PITCH AND DEPTH CONTROL 


A. OPTIMAL SOLUTION TO THE LINEARIZED SUBMARINE 
EQUATIONS 


Automatic depth and pitch control can be achieved by 
Various methods. For example, because of the fairwater 
planes naving minimal effect on pitch a freguently used 
combination is to control the depth with the fairwater 
planes and pitch with the stern planes. Another way 1s to 
set a fixed ratio between stern planes angle and fairwater 
planes angle. In this case pitch is not "controlled" but is 
used to achieve ordered depth and the planes are used to 
Change pitch. The approach selected in this thesis was to 
let the planes individuaily achieve what ever angle was 
required to meet an optimal control based on a cost function 


to be described later. 


Because or the nonlinearity of the equations of motion 
a linearized model was required to solve the equations 
reguired for an optimal control system. The required gains 
were then determined uSing the linearized equations. 
Linearizing the equations about some operating point was 
found to be impracticai because the operating point can aot 
be constant and the nonlinearities in the equations present 
a formidable problem in taking derivatives. Linearization 
WaS accomplished by dropping all nonlinear terms with the 
following justifications. Referring to appendix A all terms 
involving iwi are dropped because the linearized nodel is in 
trin. N'=1 because Uc=U. U is a meaSured guantity and tne 
eebyecerms OL anterest are mand Q,,.1.e. the problem is 
reduced to two degrees of freedon. The last result 1s 
Obtained because only the depth and pitch are parameters to 
be controiled. Simple truncation of the equations of motion 


wot Gliminate the other nonlinear terms. The linearized 
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equations of motion are as follows 


m-Zw -lig /([wW UZw/l UzZq W 
ae 
Mw/l Ty-tq) (0 | | UOMwW/i UNg/1 | | 2 
2 2 
U Zds/l U Zab Ds 
+ 
2 2 2 


U Mas/l U Mdb/1| | Db 


Mw/l Iy-Mg 








; * 2 2 
bi UZK/1l UZq W Uezas 7 U navi | | os! 
=Inv + Inv 
Gi ~~ 
; 2 2 2 2 2 | 
>| uUMw/L UNg/L | 2 U Mds/l U udp/l | Db| 
. | E= a 


The values of the inverted matrix are constants and 
when solved for can be directly substituted as constants. 
The equations were then written in a more generai form that 
would lend itself to coding for a digital computer prograa 
aS follows. 


YT A A W B 3 Ds 
12 2t iil at 
= + 
fy) A A Q B 3 Db 
12 2a lea 2e 
Lk I L. a} kL = 


a2 





let 
Determ=({ (Iy-4q) (m~Zw) -Mw Zq j 


then 
At (Lyiig) aw+ag Mw,JU/lDetero 


2 
A pL Medwt (m-2w) Aw JU/1 Detern 
AL b (ly-Mg) 2qt2q Mg JU/Determn 


Ao t MwZat (mow) ty jU/idbetera 


2 
Bh Tl (fy-ag) Zdstaqnds WU /idetern 


2 2 
Bot iwzdst (a~2w) dds JU /i Determ 


2 
Boye (yg) adbt2q MdpjJU /lDdDetern 


Ze 
By te Mwadb +(m-Zw) MdpjJU /l Detern 


The problem was treated as a linear tracking proolen 
With the control designed around the linearized mode]. The 
methods of ref. 6 were used for an optimai control with tne 


state equaticns in the general fora 


2 b<e 


= AX + BU 


wn 


which represents 


a3 





=A +B 
Q Q Db 
W 
A 
x = 
Q 
and 
DS 
ra 
U = 
Db 
A and B are asS previously defined. 


The cost function to be minimized is 


os T T 
J=1/2 J,(& OE + U RUjdt 


ns OI ON 


In the derivation of the optinai control equation the 

controi returas the system to the origin. Pf sees 

controller rc was not restricted to the origin so X was 
ot nv 


eguivaient to E where 


E= X - C 


and cr is the ordered depth and pitch. A,3,9,R, and cr are all 
constant, X and U are time varying and the plant is assumed 


to be completely controtlable. (The last assumption was 
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supported by the solution results.) T_ is considered fixed 
I 


and gate is free. The Hamiltonian is given by 
~~ 


T at iy T 
Ham = 1/2fE QE + U RU} +p AX + p BU 
~ Vv MeN fo naw 


I NIN ne AI 


The costate equation is 


~¢# Be se 
p =dHam/dX =~QX-A p + QE 


x TT 
=-QE~A p™ 


ww Pe oe Tw 


The asterisks indicate optimal controls and trajectcries. 


The Optimal control is then found by 


* ee 
dJHam/dU = RU + B p = 0 


e 
— — 


A '! ~BR B E 
= ete 2 -| |-« - 
p gor A p* 
From ref. 5 the boundary conditions are pf ea = 0 
af 


defining p =KE 


PN ee od 


taking the derivative of p and dropping the stars 


awd 


=KE+KX 


Ne AS WT 


Gs 


12) 





by substitution 


- - ih 
p =KE + KX = ~QE-A p 
~w NN rer WSIS ws fp 
- -1 7 
X = AE~BR' Bp 
“ NW AND ~/ 
pe= KE 
VJ wr 
then 
-1 7 T 
PRA + KBK BK- Q-~ ALB 
ad anv vn ~ Ww w~ wi we 
and 
-1 7 -1 T 
U =-R Bpe#=-R B KE 
a7 1/7 (ft n/ ne WI 
As described in ref. 5 Bice) = 0.0. K is Symmetric i.e. 
oe etc. so n(nt1)/2 differential equations must be 


solved where n=4. The solution of the differential equation 


for K will provide the optimal gains for JU. 


w/ 


Referring to the original equations and putting them in 


State vacriaole form 
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x 0 A 0 A x E B Ds 
2 11 21 2 11 21 
. at 0 9 1 x 0 o | |pb 
3 3 
xX 0 ry 0) A Xx Bb B 
4 | 12 22 4 12 2 
0 0 0 
A =| 0 A 0 A 
nw 11 21 
0 Q 1 
ty 12 ’ +22 
0 0 
B =!B B 
~ 11 zi 
0) 
we 22 | 
Det 0 0 0 
Q =1|0 Q 0 0 
0 0 E (an 
Oo oOo oO fo" 
ic 0 
R= 
n/ 
0 Cc 
R and Q are the weighting matrices. 


and 


The following equations result from the matrix 


manipulation. 
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11 atie 12 14 
mee=6B K +8 K 
2i ) i ar lz 2¢ 
F =B K «8B K 
31 1123 i2 34 
m= 58 K +B K 
41 1i 24 12 #4 
F =B K +B K 
le 21 12 22164 
F = K + B K 
22 2ie2 2224 
F327 =B K +B K 
23 <ie23 ze 34 
Bie =p K +B K 
2s 21 26 22 04 
° 2 2 

= = + : + 
a e Joe oe Ms 
feee= -(K +K A +X A ) + (F FP *#+F FDS 
12 11 12 ii 1¢ i2 2a iti 22 tee 
K = (F F +F F )/C 
i3 3i 42 ge ie 
feee= -(K A +K +K A ) + (F F + F FO) 
i4 lo 12 i3 14 22 41 11 42 12 

hs 

- i 2 2 
Paeee= -2(K *K A +K A) + : 
22 ( 12 ce il cé 12! an ae 17 
pee - (K Ut+thUKOUAlUlCUtTCKCOA CY) + UCFOUFOUtUFUFPOSO 
23 13 23 411 34 12 31 21 32 22 
Mee - (Ko + KK A +K.A +) - (K A +K +K_,A 3} 
24 14 z4 ii 44 i2 ez 21 23 24% 22 
F F F c 
( 41 21 42 eA 
° 2 2 
K = +P =D 
33 oer, 32 i 
Mee (K A +kK +k A ) + (PF FP +4+F F )/C 
34 23 21 33 34@ 22 > eee ae 42 32 
e 2 2 
Mme=-2(K A +t+tK..+#K A ) + (€F + F y/c 
44 2¢ 21 J3¢ 44 22 41 42 


Since only tne value of oa is known these equations must 
Fad 
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be solved backwards i.e. from t_ to on ZO.dO thrys let 
£ 
—t -t 
f£ 


ar= -dt 
when Se Po 


UAe 


= dK/dt = -dkKsarT 


“A 


The signs of all the eguations are reversed and the initial 


conditions are those conditions at eae Up to this point the 


solution 1s general and when solved the solution will 


provide gains for 


Sl | We ? 
U = -R' 3B KE 
“~ ty lv A ad 
F Je ie ef 
i121 2i 3i 41 
=-4/C E 
F F My Fr 
ie Poder 32 42 


The calculations of the gains that follecw appiy only to 
the submarine in ref. 1. The submarine described in ref. 1 
is fictitious and used for demonstration purposes only. DSL 
again is used to solve for the gains. The program used is 
program #2 which also contains the controliers designed in 


the following sections. 


Program #2 first calculates the steady state gains used 
in the ccntrolier. While the gains are calculated the ship 
Simulation is held at zero. After a specified time the 
Simulation starts with the gains determined by the specified 


weighting matrix. 


Selection or the weighting matrix to be used can be 





Started off with some judgment as to the Severity of 
different errors. For instance the pitch eyuations are in 
radians and the depth equations are in feet. It 1S obvious 
that a one radian pitch error iS more severe than a one foot 
depth error. This immediately would lead to a higher weight 
On the pitch error signal. It was found that if too nuch 
emphasis was placed on pitch tne deptn error in steady state 
would become unacceptable. This lead to a trial and error 
process to select an acceptable weighting matrix combination 


that would give the desired results. 


To operate aS an optimal controller the gains would 
have to be generated on line. This would require a 
Significant computer capability. In an infinite duration 
process the gains would be constant as demonstrated in fig. 
mrt through 7. The optimal control can then be 
approximated with a Steady state gain if the control 
interval is long with respect to the decay time of the 
Optimal gains. Pagite  fiiz |e throug wei Efi-) (indicate va 
relativeiy short rise and fal} time for the gains just 


before t . 
ife 


As determined by observing the equations for kK _ the 
gains are functions of the velocity JU. A number of 
different runs were made at different velocities and tne 


following approximate proportionalities were observed. K of 


/.. 7 
<a — eke inwersely sproportione!l to U. K  ,K_, 

~ 23 ar ; z 
3 


and K are inversely proportional to U . The runs were 


made at 25.33 feet per second so the following adjustment 


was made to the gains used in the programs. 


4O 





2 2 
2M1=(25.33) K /U 


2 2 
ZM2=(25.33) K  /U 


3 3 
2M3=(25.33) K_ o/U 


2 2 
ZM4= (25.33) K /2 


3 3 
BB5=(25.33) Kk, /u 


2. 2 
ZM6=(25.33) K >) /U 


3 3 
ZW7= (25.33) Ko /U 


The new gains are given by Z2M1 through ZM7. The B matrix 


2 
contains U so a cancellation of U takes place in the final 


eguation. 


B. AUTOMATIC PITCH AND DEPTH CONTROL 


Figure III-8 and III-9 were Simulations with C=1.0, 
D=100.0 and E=1.0 or the control input was weighted to 1.0, 
the pitch error was weighted to 100.0, and depth error 
weighted to 1.0. The ordered depth change was 10.0 feet. The 
controller initiated the proper planes but because the pitch 
was so large when ordered depth was achieved the overshoot 
pushed the control system into instability. Also the ship 
patches in the wrong direction. The initial plane response 
is to a depth error and both planes were positioned as seen 
in fig. III-10. The ship pitched up with an angle too large 


to overcome. 


The next run was made after weighting the pitch and the 


uy 





control input more heavily. The weighting watrix used was 
C=10.0, D=300.0, and E=1.0. Figure III-11 is the resultant 
depth change for a 10 foot step input in ordered depth. The 
resultant steady state depth error was .29 feet which was 


acceptabie. AS before the ship initially pitched up to 


fe) 
about 12 and then returned to a steady state angle of 


) 
+1.28 as seen in fig. f{I11-12. Figure [11-13 shows the 
fairwater and stern plane angles used to achieve the new 


° 0 
ordered depth. The maximum ~80 stern plane angle and ~-50 


fairwater plane angle are unacnievable on a real ship anda 
the rates of change ot aepth and pitch were unnacceptable. 
When a limit constraint was piaced on the plane angles the 
system went into a limit cycle. The steady state error in 
pitch was large enough to prevent the trina controller in the 


next section from detecting a ballast change to any degree 


) 
Game accuracy. Lastly the pitch variation of 12 was 


considered unacceptable for only a 10 foot depth change. 


Larger emphasis was reguired on pitch but this caused 
mie) plane angles to become larger and if the weighting on 
plane angies was increased the system became insensitive to 
either depth error or pitch error depending on the weight. 
The problem appeared to ove a result of the initial depth 
eCrrcor. To limit this effect with a large step input the 
depth and pitch error were limited aS Shown in fig. I1I-14. 
Also a larger emphasis waS required on pitch error in the 
weighting matrix. The weighting matrix chosen was C=10.0, 
D=3000.0, and E=1.0. 


In fig. ILI-13 a significant noise level is seen in 
the plane fositions. To reduce this noise the plane angle 


ordered was filtered with a Simple low pass filter. 


The tests were run again and this time a 20 foot depth 





Change was ordered. The magnitude ot the saturation limit 
can be adjusted for diiferent results . The values selected 


O 
for saturation were + two feet for depth error and + 10 for 


pitch error. Major differences to be noted in the runs are 
the reducea rate of depth change. A 20 foot depth change 


0 
took about 45 seconds. The maximum pitch assumed was ~-.8 


- fhe steady state depth error was .2 feet and the steady 


) 
State pitch error was .14 . This run was plotted in fig 
mrr—15, 16 and 17. 


The controller was tested for its ability to maintain 
depth. fo check this the ordered depth was put at zero. 


The controller maintained a steady state depth error of .16 


0 
feet and a steady state pitch error of .14 as seen in fig. 


IfI-138 and JiiI-19. The plane angles are presented in fig. 


III-20. Note that a steaay state falrwater plane angle of 


0 
-4.0 Was reguired to maintain depth. What this implies is 


that to run the simulated submarine at a zero pitch and on 
the ordered depth with zero plane angles, ballast would have 
to be added. 


The depth controiler was also reyguired to maintain 


0 
depth in aturn. A simulated rudder angle of 10 was 


added. The rudder motion is not restrained in rate for this 
Simulation. After 50 seconds the rudder angle was put to 
zero. Figure III-21 plots the roll developed as a result of 
the manipulation. The transient in pitch and depth was small 
femendicated in f19. If1-22 and 23. In particular fig. 
III-24 shows a positive fairwater and stern plane angle 


required as the submarine settles in the stern. 


In early runs it waS apparent that the submarine was 


tlight" and ballast should be added. This was based on a 
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) 
Steady state fairwater plane angle of about -4.0 with a 


near zero depth and pitch error. iimecnentol lowing run 
ballast was added to the auxiliary tank at about 8000 pounds 
per minute for 2 minutes while the planes were in automatic 
control. Figure I1I+25 and 26 show a steady state depth 


0 
error and pitch error of about ~+-.02 feet and .04 as 


re) 
compared to .16 feet and .14 with out the ballast. Also 


the steady state fairwater plane angle, Fig. IIfI-27 , shows 


Q 0 
a steady state angle of about -1.5 compared to -4.0. 


This confirms that ballast should be added to run the Ship 


in "trim", 


The depth controlier was operating satisfactorily with 
a controller that is somewhat iess than optimal for large 
transients and near optimal for small transients. The last 
problem to solve was to determine ballast error in the 


submarine and correct it. 
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IV. AUTOMATIC TRIM CONTROL 


A. METHODS OF TRIM CONTROL 


A continuous monitoring and control of the submarine 
trim 1s required. A submarine that is "heavy" runs with the 
fairwater planes on a constant positive angle. If the 
Submarine 1S too heavy the fairwater planes angle reaches 
Saturation and a large positive pitch is required to 
Maintain depth. The speed of the submarine along with the 
angle of inclination of the planes determine the lifting 
force developed by the planes. At a high speed a submarine 
that is "heavy" may be well within the control limits of the 
planes. When the speed is decreased the decreased lifting 
forces would cause larger plane angles to be ordered until a 
possible saturation of the planes could occur. A submarine 
at any speed with a constant steady state fairwater or stern 
planes angle iS experiencing a drag that is reducing the 
efficiency of the propulsion plant. These conditions must be 


corrected to efficiently operate the submarine. 


Two options for automatic trim control were considered. 
One was to determine the steady state lifting forces caused 
by excessive ballast and display the magnitude of the 
ballast error allowing a manual correction. The next was an 
automatic controller in the loop pumping or moving ballast 


as needed to maintain trin. 
The first option would be accomplished by solving the 


linearized eguations developed in previous sections in the 


forn 


Force=AW + BGO + CDs + EDb 
Moment=AW + BQ + CDs + EDD 


Ue 





Where in the steady state, if the plane angles are non zero 
the forces and moments are equivalent to the weight of 


required ballast to be moved. This option was not explored. 


The second option was accomplished using a model 
following technique. The velocities of a linearized model 
running in a computer and those of the submarine are 
compared as in fig. IV-1. Both the model and the plant have 
the same plane angle inputs. For zero plane angles and 
operation on ordered depth the velocities in the vertical 
-direction and around the Y axis for both the plant and model 
are zero. Perturbations caused by the planes to achieve a 
new ordered depth should cause both the plant and the model 
to accelerate in the same direction. However if a steady 
State plane angle 1S required to maintain depth or pitch in 
the submarine these plane angles represent accelerations in 
the model. A comparison of the submarine which is at nearly 
zero velocity and the model which is accelerating would 
indicate a ballast differential between a model in "“trin® 
and a submarine out of "trim". There will be some error in 
the model but for small perturbations the error should be 
small. The model actually used was the sane as used in 
previous sections for the derivation of the depth 


controller. 


The need for the watch officers to maintain absolute 
control of the ballast shirted in or out of the ship would 
prevent an automatic trim controller from working in other 
than a paSSive override mode. In this mode the controller 
Signals the operator that corrective action must be taken to 
Maintain the proper trim. Themeoperaton, iat conditions 
permit will actively allow the controller to perform its 
task. With no action on the part of the operator the system 


would not operate. 
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B. MODEL COMPARISON FOR TRIM CONTROL 


In all the Simulations that follow the velocities used 
are Wand Q. The error resulting is small as long aS pitch 
and roll are small. It wWaS observed that with automatic 
control the fitch was kept relatively small. The effect of 
the roll will be discussed later. All runs are made with 


the submarine in an automatic depth control node. 


To get an idea of the magnitude of the velocities 
invoived the first runs conSisted of comparing the model and 
Submarine velocities in the W and Q direction while manuaily 
Oietiating a Shift of ballast. In fig. IV-2 ballast was 
pumped into the auxiliary tank. The bottom curve is the 
fairwater planes angle which moved towards a ze Lo position. 
This is the same result observed in the previous section. 
The top curve is the stern planes angle and the niddle curve 
1s 10 times the velocity error. The error signals are 
multiplied to fit them on the same plot as the planes. Note 
that the velocity error crossed zero at approximately the 
Same time as the fairwater planeS were at zero. This was 
the desired response to adjust trim for a zero fairwater 


plane angle. 


This simulation run was repeated, this time ballast was 
shifted from after trim tank to forward trim tank. In fig. 
IV-3 the fairwater planes angle and stern planes angle, 
curve 1 and 2, were adjusted to overcome the turning moment 
caused by the ballast shift. Curve 3 is 2000 times the 
velocity errcr in the Q direction. This curve crossed zero 
at about the same time as the stern plane angle crossed 
zero. This was the desired response for a zero stern plane 


angle when the supmarine was in "trim", 


The controller used was an on-off controller. The pump 


Operated at a fixed rate. The error detectors were to sense 
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an error, determine the direction of the error and then 
Operate the pump. To prevent the pump from operating 
continuously a deadzone was inserted in the error. signal 
Channel. the deadzone was used to limit cycling of the pump 
during transients but its magnitude was constrained by the 
plane angles aliowed in steady state. Adjustment of tae 
deadzone would also determine the frequency of pump cycling 


and the steady state accuracy. 


In fig. IV-4 the vertical velocity error and the steady 


state accuracy were plotted with a history of maneuvers as 


follows. At time zero the rudder was put at a left 10. AP ake 
50 seconds it was moved back to zero. At 100 seconds the 
depth ordered was ~20.0 feet. The peak observed at time 70 
seconds would prevent any fine control of trim. This peak 
was the result of the settling of the submarine in a turn as 
Observed in ,»revious sections. The settling or squatting is 
compensated for with plane angle and these piane angles 
accelerate the model. This is not a Serious poblem and can 
be solved ina number of ways. One way would be to increase 
the degrees of freedom in the model. This would of course 
be a more ditficuit model to construct in a computer. 
Alternately the system can disengaged ina turn. This would 
reduce the capability of the system by a small degree. 
Since the time Spent in actually turning is smail 
disengaging might pe tne best course of action. For use in 
this thesis it was assumed that this second course of action 
was taken. The deadzones were then selected at .0008 
radians/second £Or rotational velocity error and .08 


feet/fsecond fcr verical veiocity error. 


The controller had two separate modes. One was to move 
ballast in and out of auxiliary tank to sea. The other was 
to move pallast between forward and after trim tanks. These 


modes were first checked out separately. First the 
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ballast the pump was shifted to pump from auxiliary for 79 
seconds. At this point the controller stopped pumping out 
of the auxiliary tank and again pumped from after trim to 
forward trim for 20 seconds and then the pump stopped. 70 
seconds after the pump stopped it started and pumped for 20 
seconds out of the auxiliary tank. The end result on the 


0 0 
Steady state angle of the fairwater planes was .6 and .2 


On the stern planes. Figure IV-9 shows the tank levels as 
ballast is moved to achieve a neutral trim. Figure IV-10 
and 11 plot the planes angles. Tne planes move toward a 


zero steady state angle and a neutral trim is approached. 


To denomstrate the effect of the deadzone the vertical 
velocity error deadzone was reduce to .04 feet/fsecond fron 
-08 feet/seccnd. Figure IV-12 is the tank levels where the 
initial out of trim condition was the same as for fig. IV-9. 
BeecoOmparison of £lg. iIV-10 and 11 with fig. IV-13 and 14 
shows the steady state planes angles were closer to zero 
With the smaller deadzone. Figure IV-15 and 16 show the 
depth controller while ballast is shifted. Pi tcho es 


0 
controlled to .02 error and depth to a .1 foot error. 


When the ship was ordered to change depth by 20 feet 
With the Ship in trim the nodel followed ciosely enough so 


that no pump commands were given. 


An additional feature in the system was tested. On 
some maneuvers it is desirable to operate the submarine 
“Nheavy". This is normally accomplished by simply adding for 
example 5000 lbs to the auxiliary tank. Whetner or not the 
Submarine is actuaily 5000 lbs heavy depends on the initial 
state of the trim. If the submarine is initially 2500 lbs 
light it is now only 2500 lbs heavy vice the desired 5000 
lbs. By adding a weight forcing term to the model in the W 


eguation and assigning a value to this forcing term when the 
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heavy or light conditon 1s desired the pump will shift water 
into or out of the auxiliary tank until the desired 
condition is reached. The accuracy of accomplishing this is 
determined by the deadzone in the vertical velocity error 
Signal. For the submarine Simulated a value of .02 


feet/second worked nicely. 


The program used for the Simulation is program #2. Once 
the gains and cofactors are determined they can be read into 
the program aS parameters and the subroutines and Riccati 
Equations (gain eyuations on page 38) can be removed. The 


controllers with filters as used is shown in fig. IV-17. 
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Figure IV-l 


Model Comparison Block Diagram 
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Figure IV-2 
Curve 1, Fairwater Planes Angle vs. Time 
Curve 2, Stern Planes Angle vs. Time 
Curve 3, Vertical Velocity Error vs. Time 
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Figure IV-3 
Curve 1, Fairwater Planes Angle vs. Time 
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Figure IV-4 
Vertical Velocity Error vs. Time 
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Figure IV-5 
Auxiliary Tank Level vs. Time 


8 3 





(.2 Feet/sec.)/inch 


2 Degrees/inch 


8.00 


6.00 


4.00 


2.00 


| Ay rat Ly yy 


A tlh Ea 
_™ NTE 


is fs’ 
al 
. 


00 


= od 
a id 


- 00 


-0. 00 6. 00 10. 00 15.00 20.00 25. 00 
(xio )* 
50 Seconds/inch 


Figure IV-6 
Curve 1, Fairwater Planes Angle vs. Time 
Curve 2, Stern Planes Angle vs. Time 
Curve 3, Vertical Velocity Error vs. Time 
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Figure IV-7 
Ferward. Jrim Tank vs. Time 
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Figure IV-8 
Curve 1, Fairwater Planes Angle vs. Time 
Curve 2, Stern Planes Angle vs. Time 
Curve 3, Rotation Velocity Error vs. Time 
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Figure IV-9 
Curve 1, Auxiliary Tank Level vs. Time 
SUivesa.  Oryward Trim Tank Level vs. Time 
Curve 3, After Trim Tank Level vs. Time 
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Figure IV-10 


Fairwater Planes Angle vs. Time 
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Figure IV-11 
Stern Planes Angle vs. Time 
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Figure IV-13 
Fairwater Planes Angle vs. Time 
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Figure IV-14 
Stern Planes Angle vs. Time 
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Figure IV-16 
Depth vs. Time 
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Controller Block Diagram 
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Ve. LISCUSSION, CONCLUSIONS AND RECOMMENDATIONS 


7 = = —_— ee = — SS SP a a = _ _ a ae ee ee ae ae SS —_— SS 


me )|)|«||CU DISCUSSION 


In the design of optimal controllers a major difficulty 
is in determining a weighting matrix that produces the 
desired response. With the use of computers many Simulations 
can be made and the weights adjusted to obtain the response 
needed. In this thesis no attempt was made to weight the 
velocities with the errors. This is an option that could 
give more control over the output response, but this would 
also greatly increase the difficulty in selecting the 


weighting matrices. 


The optimai controller designed gave excellent results 
and produced steady state errors that were smaller than can 
be normally meaSured. The magnitude of the error achieved 
is a function of the weighting matrix selected and the 


forcing inputs to the systen. 


The depth controller was not deSigned as a depth 
changer, but as a regulator. However it did function for 
depth changing inputs. The depth changing rate was 
determined by a Saturation limit in the depth and pitch 
error inputs . The limiter is reguired because the flanes 
have a finite capability and the use of optimal control 
assumes an infinite control input is available. The 
Saturation limit can be adjusted for different plane angle 
responses to large depth changes. This would change the 
rate of depth change and the pitch angle assumed during the 
transient. The magnitude of the saturation limit chosen was 
based on limiting the planes angles when a new ordered depth 
Was given and further study of any range of limits was not 


made. 
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Trim control was accomplished using a model comparison 
technique. The results of this technigue were in some 
respects better than expected. The trim controller was 
accurate in the final results and its corrective action was 
achieved earlier than a manual operator could accomplish 
because the error was detected sooner. The model was run in 
two degrees of freedom. If the model is left at two degrees 
of freedom scmething must be done to compensate for tie 
effects observed ina roll. The system would not reliably 
detect a change in trim in a turn and its ability to detect 
a change in trim during a large vertical transient is 
questionable. However, the controller would not mistake a 
vertical transient as a demand or need for ballast change. 
Because of this last result it was not necessary to 
compensate or change the system in a dynamic vertical 
Situation. The trim controller was not based on the depth 
controller used so it would conceivably work with any depth 
controller as long as the planes were used to compensate for 
pitch and depth errors. Filters would be reguired to give 


continuous average plane angles. 


Using the trim controller to adjust the submarine to a 
"heavy" or "light" condition waS not part of its intended 
function but the controller was able to accomplish this with 
reasonable accuracy. The accuracy was largely determined by 
the dead zone in the Signal channel. The deadzone would have 
to be large enough not to mistake noise or ordered depth 
transients for a trim error. This contraint imposes a limit 
on the accuracy that could be achieved. The priority 
assigned in adjusting trim was arbitrariiy assigned to first 
correct an Overall weight problem then adjust for a fore and 


aft trim. 
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B. CONCLUSIONS 


Optimal control techniques can be used to design a 
depth controller that achieves a steady state condition as 
accurately as desired. Adjustment of the weighting matrices 
would emphasize those errors that are most important and 
also adjust the amount of response py the planes to 
accomplish the control. The insertion of nonlinearities (in 
this case a Saturation limit) would allow the controller to 
give a reasonable response to any input and provide for a 
depth changing capability as well as the depth maintaining 


Capability designed for. 


Steady state gains were used instead of time varying 
gains, and this did not degrade the functioning of the 
controller but did provide a significant simplification in 
the controiler. A controller designed with the metnods 
described would provide a vast improvement over manual 
Operation. The gains provided are only a function of speed 
which makes realization of the controller relatively simple. 
The contrcller output itself is also a function of pitch 
rate and depth rate. This provides an anticipatory feature 


to the contrcller. 


In the past trim control has been purposely left to 
manual control. This was for safety reasons and because the 
detection of an out of trim condition was a matter of 
judgement based on the observation of plane positions and 
theretore only good in the steady state. The passive 
override mode discussed in previous sections removes the 
safety constraints in the use of an automatic trim 
controlier. The simulations run demonstrate that detection 
of trim error can be accomplished accurately in a dynamic 
Situation. The method still uses the planes positions for 
its determination but the lifting forces of both sets of 


planes and the turning moments of both sets of planes are 
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used in the determination of trim. Traditionally only the 
lifting forces of the fairwater planes and the turning 
moment of the stern planes were considered for the trin 
error estimate. The controller made accurate estimates as 
seen in fig. IV-9 and 12 where it was noted that the 
controller never had to go back and correct an adjustment 


once it was made. 


The additionai feature in the trim controller of being 
able to delikerately order an out of trim condition provides 
for more poSitive control in the trim of the submarine. 
This control 1s a function of the actual submarine 
conditions and not of the assumed conditions. if athe 
conditions change during the submarine operation such as an 
ascent the trim controller wouid automatically maintain the 
out of trim condition ordered. A difficult feat for a 


Manual operator. 


C. RECOMNENDATIONS FOR FURTHER WORK 


There are six specific areas that need further study 


and work. They are enumerated in the following paragraphs. 


1. The selection of the saturation limits in the error 
Channel was arbitrary and the only criteria was to prevent 
the planes angles from becoming excessive in a depth 
changing situation. As in the actual controller there must 
be some optimum limit that would meet a specific criteria in 
changing depth. this limit could be found using a function 
minimization program where numerous Simulation runs are nade 
and compared to determine a cost based on meeting the 
criteria specified. A minimum cost would define a match to 


the criteria specified. 


2. To achieve a zero steady state error and at the sane 
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time make the controller an optimal depth changing 
controller a method of decoupling is needed. The methods of 
steady state state decoupling through the use of cascade 
compensators described in ref. 7 might provide such a 
controller with a function minimization procedure to meet a 
Specified criteria. This approach would allow for zero 


Steady state error in pitch and depth under most conditions. 


3. the controllers in this thesis were designed for 
continuous’ time. To provide for use in a time shared 
computer environment a discrete version of the controller is 
required. For the trim controller this would mean a 
discrete modei used for comparison with the submarine and a 


Sampling of the submarine parameters. 


4. The control inputs considered for the depth 
controller were the fairwater and stern planes. The control 
inputs could be expanded to include a vertical forcing 
function and a rotational forcing function around the Y 
axis. These inputs wouid of course be the shifting of 
ballast in the trim tanks. in this case the planes angles 
would be included as states. The purpose of this scheme is 
to control the planes to a zero angle by the shifting of 
patiast. 


5. The trim controller does not work inaturn. A 
Suggesticn to correct this would be to increase the degrees 
of freedom of the model. The model would have to be 


developed ana tested. 


6. Finally, the options allowed for the pump were very 
limited and for example if the submarine is "heavy" forward 
and at the same time “heavy" overall the designed controller 
would pump ballast from auxiliary tank then pump ballast 
from forward to aft. A more efficient maneuver would be to 


pump ballast form forward trim to sea. Even for only three 
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variable ballast tanks the pumping operation would become 


Significantly more complicated. 
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EQUATIONS OF MOTION 


The following set of equations are referred to a bodv 
fixed system of axes which are coincident with the prin- 
cipal axes of {nertia of the body. The origin of this 
axis-system is located at the assumed center of mass of 
the bodv 


Equation of Motion Along the Body Axis System x-Axis 
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Equation of Motion Along the Body Axts System y-Axts 
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Equation of Motion Along the Body Axis System z-Axis 
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Equation of Motion About che Body Axis System x-Axis 
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Equation of Motion About the Body Axis System y-Axis 
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E;Squation of Motion About the Body Axis System z-Axis 
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+ 24° N Pee wvit 


+ UW. X,; COs 6 sind 


fi Multiply by ~ for large 
angles of ial near -90° 
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AUXILIARY EQUATIONS 


pt i sin8 


&: 
i) 


s) = (q - Ycos@ sind )/ cos 


y = ele ON con 0 -oa 
+w (sin sind + cos® sin@ cosh) 


+w (cos sin®@ sin» - sin® cos) 


Zo = -usinO+vcos6 sind + wcecos8 cos 


U = fas ee ew? 


Wedae By 1] al tas nt + ay'n'? | 


Xo = ucos cosh + v (sind sin@ cos - cos® sing) 


Yo = ucos@ sin~+v(cosOcos + sind sin®@ sind) 


2 when k, < n' 
= Bru] bi! + bint + byt n'? | when k <n'< k 
= Brtu?[ att cgint + coin’? | when kj <n'<ck 
- See ak + d,'n' + dy'n'? | when n'! < ky 
’ ’ : 
Gy), 5, coc) 
i ey b3 Sets of non-dimensional coeffictents used in the pro- 


pulsion equation above. ‘The set which will be {in erfect 


Cl, Co, CQ et any tite duriue a similated maneuver will depend on tne 


q e ‘ 
dy; do, ap 
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value of n’ and the nunbers Ky kok : 


NOMENCLATURE 


All symbols used in the equations of motion and in the 
auxifliary ecuatitons and relationships which appear in this 
report are defined below. Any dimensions involved will be 

_consistenet with the foot-ponund-second system of units. All 
angles are in degrees. The Fortran variables corresponding 


to these symbols are shown in Appendix g . 
SYMBOL DEFINITION 
e A dot over any symbol signi- 
fies differentiation with res- 


pect to time. 


B Buoyancy force which ts posi- 


tive upwards. 


m Mass of the submarine including 


the water in the free fisordiag 


spaces. 
Q Overall length of the submarine 
U Linear velocity of origin of 


body axes relative to an earth- 


fixed axis system. 


u Component of U along the body 
x-axis. 
Vv Component of U along the body 
y-axis. , 
e 
w Component of U along the body 
z-axis. 
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Command speed: A steady value 
of u for a given propeller rpm 
whenad, B and control surface 

angles are zero. Sign changes 


with propeller reversal. 


Longitudinal axis of the body 


fixed coordinate axis system. 


Transverse axis of the body 


fixed coordinate axis system. 


Vertical axis of the body fixed 


coordinate axis system. 


Distance along the x, axis of an 


earth-fixed axis system. 


Distance along the e axis of an 


earth-fixed axis system. 


Distance along the z axis of an 


earth-fixed axis system. 


Component of angular velocity 


about the body fixed x-axis. 


Component of angular velocity 


about the body fixed y-axis. 


Component of angular velocity 


about the body fixed z-axis. 


The z coordinate of the center 
of buoyance (CB) of the subma- 


rine. 





Sb, De 


Sr ) Dr 


55 , Ds 


3 


rm && eG O 


x 


ti 


Angic of attack. 
Angle of dritc. 


Deflection of bowplane ( or 


eailplane ) 

Deflection of rudder. 
Deflection of sternplane. 

The ratio u./u. 

Pitch angle. 

Yaw angle. 

Roll angle. 

Mass density of sea water. 
Weight of water blown from a 
particular ballast tank {dent- 
ified bY the integer lassigned 
to the index 1. 

Angular velocity. 


Time. 


Location along the body x-axis 


of the center of mass of thé pth 


ballast tank when this tank is 


filled with sea water. 
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M.’, M ' 
q rr 
f 
M , M 
Iwiq? oe 

M t 


( 


Fx) p 


EZ 


Propulsion force (see auxt- 


liary equations and relation- 


ehips). 


Monent 


marine 


Moment 


marine 


Moment 


marine 


of tnertia of a sub- 


about the x-axis. 


of tnertia of a sub- 


about the y-axis. 


of inertia of a sgub- 


about the z-axis. 


Non-dimensional constants each 


of which is assigned to a partit- 


cular force term in the equation 


of motion about the body x-axts. 


Non-dimenstonal constants each 


of which is assigned to a parti- 


cular force term in the equation 


of motion about the body y-axis. 


we at see Woo ero Mire 
a Nye Noone ea Nao 5 ae 
aeiv| ° Ge ; eat Sune : 

a 


ps Be tr tp u vr wq 
Kou? x? ue roi x aed 
XSbsb° vvn ° feat Apa canat 
SerSrn' 

“og om Pq’ "lel? ee : Yop 
*vdr! oY Yieler ? Vea as 
Bae ty Yyivl ee oat 
Dent | eee se Soe 
aa oe ie aE aur? ye 
B2yp's 2q'* Males? *wlal * 
Zon'? yo Zo al a : ate : 
“ope Deve Ze zen? ea a 
ho AGce yo SUAS 


i] 
wylvla 
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Non-cimenstfonal constants each 
of which is asaigned to a parti- 
cular force term in the equation 


of motion about the body z-axis. 


Non-dimenstfonal constants cach 
of which is assigned to a parti- 
cular force term in the equation 


of motion along the body x-axis. 


Non-dimenstonal constants each 
of which is essigned to a parti: 
cular force term in the equatlo 


of motion along the body y-axis 


Non-dimensional constants each 


e 
of which is assigned to a parti- 


cular force term in the equatio- 


of motion along the body z-axis 
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